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Abstract. Using the potential model and thermodynamical quantities obtained in 
lattice gauge calculations, we determine the spontaneous dissociation temperatures of 
color-singiet quarkonia and the 'quark drip lines' which separate the region of bound 
QQ states from the unbound region. The dissociation temperatures of J/ip and Xb in 
I quenched QCD arc found to be 1.62Tc and l.lSTc respectively, in good agreement with 

■"sj" ■ spectral fmiction analyses. The dissociation temperature of J/ip in full QCD with 2 

, flavors is found to be lA2Tc. For possible bound quarkonium states with light quarks, 

(—^ ' the characteristics of the quark drip lines severely limit the stable region close to the 

■ phase transition temperature. Bound color-singlet quarkonia with light quarks may 

I exist very near the phase transition temperature if their effective quark mass is of the 

^ ■ order of 300-400 MeV and higher. 



PACS numbers: 25.75.-q 25.75.Dw 



1. Introduction 



^ I The degree to which the constituents of a quark-gluon plasma can combine to form 
5^ \ composite entities is an important property of the plasma. It has significant implications 
on the plasma equation of state, the probability of recombination of plasma constituents 
prior to the phase transition, and the chemical yields of the observed hadrons. Recent 
spectral analyses of quarkonium correlators indicated that J/ip may be bound up to 
l.6Tc where Tc is the phase transition temperature P3I2]- Subsequently, there has been 
renewed interest in quarkonium states in quark-gluon plasma as Zahed and Shuryak 
suggested that QQ states with light quarks may be bound up to a few Tc jS] • Quarkonium 
bound states and instanton molecules in the quark-gluon plasma have been considered 
by Brown, Lee, Rho, and Shuryak [4 . As heavy quarkonia may be used as diagnostic 
tools [S], there have been many recent investigations on the stability of heavy quarkonia 

in the plasma jniiziiHiiniiiniiniiiaiisiiiniiiHiiii. 

Previously, DeTar [T7], Hansson, Lee, and Zahed [THj, and Simonov jTHj observed 
that the range of the strong interaction is not likely to change drastically across the 
phase transition and suggested the possible existence of relatively narrow low-lying QQ 
states in the plasma. On the other hand, Hatsuda and Kunihiro |20ij considered the 
persistence of soft modes in the plasma which may manifest themselves as pion-like 
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and sigma-like states. The use of the baryon-strangeness correlation and the charge 
fluctuation to study the abundance of light quarkonium states in the plasma have been 
suggested recently j2Il I22j • 

We would like to use the potential model to investigate the composite properties 
of the plasma and to determine its 'quark drip lines'. Here we follow Werner and 
Wheeler [21] and use the term 'drip line' to separate the region of bound color-singlet QQ 
states from the unbound region of spontaneous quarkonium dissociation. It should be 
emphasized that a quarkonium can be dissociated by collision with constituent particles 
to lead to the corresponding 'thermal-dissociation line' and 'particle-dissociation line', 
which can be interesting subjects for future investigations. 

2. The color-singlet Q-Q potential 

The most important physical quantity in the potential model is the Q-Q potential 
between the quark Q and the antiquark Q in a color-sinbglet state. Previous works in 
the potential model use the color-singlet free energy Fi jHl El HE] or the color-singlet 
internal energy Ui jTUl El US] obtained in lattice gauge calculations as the color-singlet 
Q-Q potential without rigorous theoretical justifications. The internal energy Ui is 
significantly deeper and spatially more extended than the free energy Fi. Treating 
the internal energy f/i as the Q-Q potential led Shuryak and Zahed to suggest the 
possibility of bound color-singlet quarkonium states with light quarks in the plasma 
|3]. The conclusions will be quite different if one uses the free energy Fi as the Q-Q 
potential. 

While Fi or f/i can both be used as the Q-Q potential at T = (at which Fi = Ui), 
the situation is not so clear in a thermalized quark-gluon plasma. In Ref. jE], we show 
that the equation of motion for the QQ state can be obtained in two steps. First, 
one considers a lattice gauge calculation for a static pair of color-singlet Q and Q at 
a fixed separation R at the temperature T and obtains the free energy Fi{R), the 
internal energy Ui{R), and their difference TS'i(r) = Ui — Fi. In the second step, one 
considers a Q and a. Q in dynamical motion in the quark-gluon plasma. The dynamical 
degrees of freedom can be taken to be the set of quarkonium and plasma constituent 
wave functions and their corresponding state occupation numbers. The equilibrium 
of such a thermalized system at a fixed temperature and volume occurs when the 
grand potential is a minimum. One explicitly writes out the grand potential, which 
is the sum of the free energy and appropriate Lagrange multiplier terms, in terms of 
the quarkonium and plasma constituent wave functions, the state occupation numbers, 
and the mutual Q-Q, (^-(constituent), Q- (constituent), and (constituent)- (constituent) 
interactions. The minimization of the grand potential with respect to a color-singlet 
quarkonium wave function leads to the Schrodinger equation for the quarkonium 
state with a color-singlet potential U^q containing only Q-Q, (5-( constituent), and Q- 
(constituent) interactions. On the other hand, the internal energy Ui in the lattice 
gauge calculation contains not only contributions from these Q-Q, (^-(constituent), and 
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constituent) interactions, but also an additional contribution from the (constituent)- 
(constituent) interaction. Therefore, the color-singlet potential U^q is determined by 
the internal energy Ui after subtracting out the quark-gluon plasma internal energy. 
These results are further supported by the presence of a similar relationship between 
the total internal energy and the Q-Q potential in the analogous case of Debye screening, 
for which analytical expressions for various thermodynamical quantities can be readily 
obtained and compared 



3. The color-singlet Q-Q potential in quark-gluon plasma 

In order to subtract out the /^-dependent internal energy of the quark-gluon plasma 
from the total internal energy to obtain the color-singlet Q-Q potential in the plasma, 
additional lattice gauge calculations may be needed to evaluate the quark-gluon plasma 
internal energy in the presence of a color-singlet Q and Q pair. It is nonetheless useful 
at this stage to suggest approximate ways to evaluate this quark-gluon plasma internal 
energy. The subtraction can be carried out by noting that the quark-gluon plasma 
internal energy density e is related to its pressure p and entropy density a by the First 
Law of Thermodynamics, e = Ta — p, and the quark-gluon plasma pressure p is also 
related to the plasma energy density e by the equation of state p{e) that is presumed 
known by another independent lattice gauge calculation. Thus, by expressing p as 
(3p/e)(e/3) with the ratio a(T) = 3p/e given by the known equation of state, the plasma 
internal energy density e is related to the entropy density Ta by e = [3/(3 + a{T))]Ta, 
and the plasma internal energy integrated over the volume is given by [3/ (3 + a(T))]TS'i 
where TSi has already been obtained as Ui — -Fi. The proper Q-Q potential, U^q^ 
as determined from Ui by subtracting the plasma internal energy (which is equal to 
[3/(3 + a{T))]{Ui — Fi)), is then a linear combination of Fi and Ui given by 6J, 

The potential U^q is approximately Fi near Tc and 3Fi/4 + t/i/4 for T > 1.5Tc jO]. 

In the spectral function analyses, the widths of many color-singlet heavy 
quarkonium broaden suddenly at various temperatures [H 12 US] • In the most precise 
calculations for J /if) in quenched QCD using up to 128 time-like lattice slices, the 
spectrum has a sharp peak for 0.78Tc <T< 1.62Tc and a broad structure with no sharp 
peak for 1.70Tc < T < 2.33Tc [1,. The spectral peak at the bound state has the same 
structure and shape at 0.78Tc as it is at 1.62Tc. If one can infer that J/ip is stable and 
bound at 0.78rc, then it would be reasonable to infer that J/ip is also bound and stable at 
1.62Tc. The spectral function at 1.70Tc has the same structure and shape as the spectral 
function at 2.33Tc. If one can infer that J/tp is unbound at 2.33Tc, then it would be 
reasonable to infer that J/tp become already unbound at 1.70 Tc. Thus, from the shape 
of the spectral functions, the temperature at which the J/ ip width broadens suddenly 
from 1.62rc to 1.70Tc corresponds to the J/ip spontaneous dissociation temperature. 
Dissociation temperatures for Xc and Xb in spectral analyses in quenched QCD have 
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also been obtained j2l C^- We list the heavy quarkonium spontaneous dissociation 
temperatures obtained from spectral analyses in quenched QCD in Table I. They can 
be used to test the potential model of U^q, the linear combination of Ui and Fi proposed 
in Eq. (Q), as well as Fi and Ui. 



Table I. Spontaneous dissociation temperatures obtained from different potentials. 





Quenched QCD 


Full QCD (2 flavors) 


States 


Spectral 
Analysis 




Fi 


Ui 




Fi 


Ui 




1.62-1. TOT] 


1.62 


1.40 Tc 


2.60 


1.42 Tc 


1.21 Tc 


2.22 Tc 


Xc 


below I.IT^^ 


unbound 


unbound 


LIST, 


1.05 Tc 


unbound 


1.17Tc 






unbound 


unbound 


1.23 T, 


unbound 


unbound 


1.11 Tc 






4.1 T, 


3.5 T, 


~5.or. 


3.40 Tc 


2.90 Tc 


4. 18 Tc 


Xb 


1.15-1.54T| 


1.18T, 


1.10 Te 


1.73 


1.22 Tc 


1.07 Tc 


1.61 Tc 






1.38 Te 


1.19T, 


2.28 


1.18 Tc 


1.06 Tc 


1.47 Tc 



tRef.p, ^Ref.[21, »Ref.ll3| 



4. Quark drip-lines in quark-gluon plasma 

To evaluate the Q-Q potential, we use the free energy Fi and the internal energy Ui 
obtained by Kaczmarek et al. in both quenched QCD JT] and full QCD with 2 flavors 
|12j . In quenched QCD, Fi and Ui can be parametrized in terms of a screened Coulomb 
potential with parameters shown in Figs. 1 and 2 of Ref. [Hj. In full QCD with 2 flavors, 
Fi and Ui can be represented by a color-Coulomb interaction at short distances and a 
completely screened, constant, potential at large distances as given in Ref. |7], although 
other alternative representations have also been presented fSl ^] . To determine the 
U^qq potential as given by Eq. (P), we also need the ratio a(T) = 3p/e from the plasma 
equation of state. We use the quenched equation of state of Boyd et al. |2S1 for quenched 
QCD, and the equation of state of Karsch et al. [2^1 for full QCD with 2 flavors. Using 
quark masses m^. = 1.4 GeV and ruf, = 4.3 GeV, we can calculate the binding energies 
of color-singlet heavy quarkonia and their spontaneous dissociation temperatures. We 
list in Table I the heavy quarkonium spontaneous dissociation temperatures calculated 
with the Uqq potential, the Fi potential, and the Ui potential, in both quenched QCD 
and full QCD. 

The spontaneous dissociation temperatures of J/ijj and Xb obtained with the U^q 
potential in quenched QCD are found to be 1.62Tc and 1.18Tc respectively. On the other 
hand, spectral analyses in quenched QCD give the spontaneous dissociation temperature 
of 1.62-1.70Tc for J/^ P and 1.15-1.54Tc for Xb HB]. Thus, the U^^^ potential of Eq. (0) 
gives spontaneous dissociation temperatures that agree with those from spectral function 
analyses. This indicates that the Uqq potential, defined as the linear combination of 
Ui and Fi in Eq. ((T)), may be an appropriate Q-Q potential for studying the stability of 
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Figure 1. The drip lines in quenched QCD calculated with the f^gg: ^i, and 
Ui potentials. The solid-circle symbols represent results from lattice gauge spectral 
function analyses. 



5. The stabihty of a Q-Q pair 

To examine the stabihty of a Q-Q pair, we consider the quark mass uiq as a variable and 
evaluate the spontaneous dissociation temperature as a function of the reduced mass 
/ired = i^Qi^Q / ijnq + i^q)- The results for quenched QCD are shown in Fig. 1. A state 
is bound in the (T/T^, /ired) space above a drip line and is unbound below the drip line. 

The drip lines for quenched QCD in Fig. 1 are useful to provide a comparison of 
the potential model results with those from spectral analyses in the same quenched 
approximation. However, as dynamical quarks may provide additional screening, it is 
necessary to include dynamical quarks to assess their effects on the stability of quarkonia. 
Accordingly, we use the U^q potential defined by Eq. ((T)) with Fi, f/i, and a{T) evaluated 
in full QCD with 2 fiavors 122] to determine the drip lines shown as solid curves in 
Fig. 2. In comparison with quenched QCD results, the Is drip line in full QCD is shifted 
to lower temperatures for T > 1.2Tc while the Ip drip line in full QCD is only slightly 
modified. 
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Figure 2. The Is and Ip quark drip lines calculated with the C/lll potential in 2-flavor 
QCD (solid curves) and quenched QCD (dashed curves). 

For heavy quarkonia, the results in Table I and Fig. 2 indicate that for the U^q 
potential J/ip, T, and Xh may be bound in the plasma up to 1.42Tc, 3.40Tc, and 
1.22Tc respectively. As light quarks have current-quark masses of a few MeV, we may 
be advised that non-relativistic potential models should not be used to describe hght 
quarkonia. However, due to its strong interaction with other constituents, a light quark 
becomes a dressed quasiparticle and acquires a large quasiparticle mass. In the low 
temperature region when spontaneous chiral symmetry breaking occurs with {ipip) 7^ 0, 
the quasiparticle mass is ruq ~ [ | (7 ('^^/') | + (current quark mass)], where g is the strong 
coupling constant and ('i/'V') the quark condensate [201 123 12H1 ■ This quasiparticle mass 
is the origin of the constituent-quark mass in non-relativistic constituent quark models 
j2nil2Hl- In the high temperature perturbative QCD region, the quasiparticle mass is 
rUq ~ gT/\/Q^ which is of the order of a few hundred MeV |29] . 

As the restoration of chiral symmetry is a second order transition, {iptp) decreases 
gradually as the temperature increases beyond T^. The light quark quasiparticle mass 
associated with {ipip) will likewise decrease gradually from the constituent-quark mass 
to the current-quark mass when the temperature increases beyond T^.. This tendency 
for the quasiparticle mass to decrease will be counterbalanced by the opposite tendency 
for the quasiparticle 'thermal mass' to increase with increasing temperature. In the 
region of our interest, Tc < T < 2Tc, various estimates give the light quark quasiparticle 
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masses from 0.3 GeV to 1.2 GeV [201 ED E21 ESI- As in the case of T = 0, where 
hght quarks with a constituent-quark mass of about 350 MeV mimic the effects of chiral 
symmetry breaking and non-relativistic constituent quark models have been successfully 
used for light hadron spectroscopy |2Hl Elj , so the large estimated quasiparticle mass 
(from 0.3 to 1.2 GeV) may allow the use of a non-relativistic potential model as an 
effective tool to estimate the stability of light quarkonia at Tc < T < 2Tc. At this stage 
when the uncertainties in the quasiparticle mass are much greater than the effects due 
to relativistic kinematics, an estimate based on a non-relativistic model suffices. It will 
be of interest to investigate the relativistic effects in the future, when the quasiparticle 
masses are more definitively determined. 

By examining the effects of the light quark quasiparticle masses on the quark-gluon 
plasma equation of state, Levai et al. jHOj, Szabo et al. jHl], and Ivanov et al. jH2] 
estimate that ruq is about 0.3 to 0.4 GeV at Tc < T < 2Tc. From the results in Fig. 2 
for the potential, we can estimate that as a quarkonium with light quarks has a 
reduced mass of 0.15-0.2 GeV, it may be bound at temperatures below (1.05 — 1.07)Tc. 
An open heavy quarkonium with a light quark and a heavy antiquark has a reduced 
mass of about 0.3-0.4 GeV and may be bound at temperatures below (1.11 — 1.19)Tc. 

Another lattice gauge calculation gives m^/T = 3.9 ± 0.2 at 1.5Tc [33 , which 
implies that at T=1.5Tc (or about 0.3 GeV), the quark mass will be ~1.2 GeV for 
(m, d, s) quarks. Such a 'light' quark quasiparticle mass appears to be quite large and 
may be uncertain, as the plasma will have a relatively large abundance of charm quarks 
and antiquarks, and there may be difficulties in reproducing the plasma equation of 
state. With this mass, a 'light' quarkonium will have a reduced mass of 0.6 GeV, and 
the quarkonium may be bound at temperatures below ~1.31Tc. 

In either case, the drip lines of Fig. 2 calculated with the U^q potential for full 
QCD with 2 flavors do not support bound QQ states with light quarks beyond 1.5Tc. 
A recent study of baryon-strangeness correlations suggests that the quark-gluon plasma 
contains essentially no bound QQ component at 1.5Tc [21]. 

In conclusion, we have used the thermodynamical quantities obtained in lattice 
gauge calculations to determine the quark drip lines in quark-gluon plasma for full 
QCD with two flavors. The characteristics of the quark drip lines severely limit the 
region of possible quarkonium states with light quarks to temperatures close to the 
phase transition temperature. If the light quark quasiparticle mass is 0.3-0.4 GeV as 
estimated from jHOlEIlESI, a quarkonium with light quarks may be bound below (1.05- 
1.07)Tc and an open heavy quarkonium may be bound below (l.ll-1.19)Tc. 
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